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PACS 75. 10. Nr - Spin-glass and other random models 

PACS 72. 10. Fk - Scattering by point defects, dislocations, surfaces, and other imperfections (in- 
cluding Kondo effect) 
PACS 72 . 15 . Jf - Thermoelectric and thermomagnetic effects 

Abstract. - We study the electrical resistivity in the quantum Griffiths phase associated with 
the antiferromagnetic quantum phase transition in a metal. The resistivity is calculated by means 
of the semi-classical Boltzmann equation. We show that the scattering of electrons by locally 
ordered rare regions leads to a singular temperature dependence. The rare-region contribution to 
the resistivity varies as T with temperature T, where A is the usual Griffiths exponent which takes 
the value zero at the critical point and increases with distance from criticality. We find similar 
singular contributions to other transport properties such as thermal resistivity, thermopower and 
the Peltier coefficient. We also compare our results with existing experimental data and suggest 
new experiments. 
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Introduction. — Quantum phase transitions [T] occur 
at zero temperature when an external parameter such as 
magnetic field, pressure or chemical composition is varied. 
They are driven by quantum rather than thermal fluc- 
tuations. At continuous quantum phase transitions, i.e., 
quantum critical points, the quantum fluctuations driv- 
ing the transition diverge and become scale invariant in 
space and time. These fluctuations dominate the mate- 
rial's properties in the vicinity of the quantum critical 
point at low but non-zero temperatures. In metallic sys- 
tems, they can cause strong deviations from the conven- 
tional Fermi-Liquid behavior of normal metals [2] . 

Impurities, defects or other kinds of quenched disorder 
can significantly modify the low temperature behavior of 
quantum many-particle systems. The interplay between 
dynamic quantum fluctuations and static disorder fluctu- 
ations leads to much more dramatic effects at quantum 
phase transitions than at classical thermal phase transi- 
tions, including quantum Griffiths singularities [3-5|, in- 
finite randomness critical points featuring exponential in- 
stead of power-law scaling [6j[7j and the smearing of the 
phase transition |8j. These unconventional phenomena 
are caused by large spatial regions (rare regions) that 
are devoid of impurities and can show local order even 
if the bulk system is in the disordered phase. The fluc- 
tuations of these rare regions are very slow because they 
require changing the order parameter in a large volume. 
Griffiths showed that this leads to a singular free en- 



ergy in a whole parameter region which is now known 
as the Griffiths phase. The probability V{L d ) for finding 
an impurity-free rare region with linear size L in a dis- 
ordered system is exponentially small in its volume L d , 
V(L d ) ~ exp(— cL d ) with c being a constant that depends 
on the disorder strength. In systems in which the char- 
acteristic energy e of such a rare region decays exponen- 
tially with its volume, e ~ exp(— bL d ), the resulting den- 
sity of states is of power-law type, p(e) cx e , where 
A = c/b is the non- universal Griffiths exponent. It varies 
systematically within the Griffiths phase and vanishes at 
the critical point. The power- law density of states pie) 
leads to non-universal power-law quantum Griffiths singu- 
larities of several thermodynamical observables including 
the specific heat, C ~ T A , and the magnetic susceptibil- 
ity X ~ T ■ The zero-temperature magnetizatio n- fie ld 
H x 



curve behaves asM~ H (for reviews, see Refs. |9j|10j). 

Quantum Griffiths phases have been predicted to occur 
not only in localized magnets but also in metallic systems 



11 



13 , but clear-cut experimental verifications have been 



absent for a long time. Only recently, quantum Griffiths 
phases have been observed in experiment in a number of 
systems such as magnetic semiconductors [l4r[l6], Kondo 



nets 



lattice ferromagnets 17|18 and transition metal ferromag- 



19 . The lack of experimental evidence for quantum 
Griffiths phases in metals may be (at least partially) due 
to the theories being incomplete: while the thermodynam- 
ics in quantum Griffiths phases is comparatively well un- 
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derstood, very little is known about the experimentally 
important and easily accessible transport properties. 

In this Letter we therefore study the electrical resistiv- 
ity in the quantum Griffiths phase of an antiferromagnetic 
metal by means of the semi-classical Boltzmann equation 
approach. In the same manner, we also investigate other 
transport properties such as the thermal resistivity, the 
thermopower and the Peltier coefficient. We find that 
the scattering of the electrons by spin-fluctuations in the 
rare regions leads to singular temperature dependencies 
not just at the quantum critical point but in the entire 
antiferromagnetic quantum Griffiths phase. The rare re- 
gion contribution to the resistivity varies as Ap oc T x 
with temperature T, the contribution to thermal resistiv- 
ity behaves as AW oc T A_1 , and the thermopower and the 
Peltier coefficient behave as AS* oc T x+1 and All oc T A+2 , 
respectively. 

Model and method of solution. Let us now 

sketch the derivation of these results. The transport prop- 
erties of the itinerant antiferromagnetic systems we are 
interested in can be described by a two-band model con- 
sisting of s and d electrons 20 21 . The Hamiltonian has 
the form H = H s + + H s _ ( i, where H s and are 
the Hamiltonians of s and d electrons, respectively. H s -d 
corresponds to the exchange interaction between s and d 
electrons. Only the s electrons contribute to the transport 
properties. They are scattered by the spin-fluctuations of 
the d electrons which are assumed to be in the antiferro- 
magnetic quantum Griffiths phase. The contribution to 
the resistivity from the scattering by the spin-fluctuations 
stems from the s — d exchange interaction term of the 
Hamiltonian 



dr s(r) • S(r) , 



(1) 



where g is the coupling between s and d electrons, s and S 
are the spin densities of the s and d electrons, respectively. 

Close to an antiferromagnetic transition in three- 
dimensional space, transport properties can be treated 
within a semi-classical approach using the Boltzmann 
equation because quasiparticles are still (marginally) well 
defined. For simplicity, we also assume that the spin- 
fluctuations are in equilibrium, i.e., we neglect drag ef- 
fects. This approximation is valid if the system can lose 
momentum efficiently by Umklapp or impurity scattering 
as is the case in a dirty antiferromagnetic system. The lin- 
earized Boltzmann equation in the presence of an electric 
field E and a temperature gradient VT, but zero magnetic 
field can be written as (22l 



v k dT 



(2) 



where /jP is the equilibrium Fermi-Dirac distribution func- 
tion. The first and second terms correspond to the rate 
changes of the electron distribution function / k due to the 



diffusion and electric field E, respectively. The last one is 
the collision term. Let the stationary solution of the Boltz- 
mann equation be / k = f?— $ k (<3/ k /$£ k ) ; where $ k is a 
measure of the deviation of the electron distribution from 
equilibrium. Then the linearized scattering term due to 



the spin-fluctuations has the form 21 23 



\ 9t J scatt 



2_l 
T 



k' 



Im X (k - k',£ k - e k ')($k - $k') ' 
^E^k'(s k -^k')(*k- $ k') ( 3 ) 



T 



k' 



where n(e k — £ k ') is the Bose-Einstein distribution func- 
tion and x is thfi total dynamical susceptibility of the spin- 
fluctuations of the d electrons. 

Electrical resistivity. — In order to calculate the 
electrical resistivity we consider Ziman's variational prin- 
ciple 22 1 . The resistivity p is given as the minimum of a 
functional of $ k [22] p] 



P[* k ] 



i //(* k -* k ') 2 rE'dkdk' 



2T 



Cf v k $ 



^dk) 



(4) 



where 



dw Pk>(oj)5(sk< — £k + ui) . 



(5) 



Quantum Griffiths effects in disordered metallic systems 



are realized both in Heisenberg magnets 13 14 and in 



Ising magnets. In the latter case, they occur in a transient 
temperature range where the damping is unimportant |11| . 
In the following, we consider both cases. 

As we are interested in the rare-region contribution to 
the resistivity in the Griffiths phase, we need to find the 
rare region dynamical susceptibility which is simply the 
sum over the susceptibilities of the individual rare regions. 
The imaginary part of the dynamical susceptibility of a 
single cluster (rare region) of characteristic energy e in the 
quantum Griffiths phase of a disordered itinerant quantum 
Heisenberg antiferromagnet is given by 



Imxd(q,w;e) = 



2 (T) + 7 2 cj 2 



F?(q), 



(6) 



where p is the moment of the cluster and 7 is the damping 
coefficient which results from the coupling of the spin- 
fluctuations and the electrons. e(T) plays the role of 
the local distance from criticality. For high temperatures 
7T > e, e(T) fa T and for low temperatures 7T < e, 
e(T) ps e. F e (q) is the form factor of the cluster which 



x We set Plank's constant, electron's charge and Boltzmann con- 
stant h = e = kg = 1 in what follows. 
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encodes the spatial magnetization profile. For random 
quantum Ising models the imaginary part of the dynami- 
cal magnetic susceptibility of a single cluster (rare region) 
is given by 



Imxd(q,w;e) 



7T— tanh | — ) 
4 \2TJ 



x [S(e-u)-S(e + uj)]F 2 (ci). (7) 



To get the total rare-region susceptibility, we integrate 
over all rare regions using the density of states p(e), 



Imx(q,w) 



dep(e)Imxd(q,w;e) 



(8) 



where A is an energy cut-off. The precise functional form 
of F e (q) is not known, since every cluster has a different 
shape and size. However, we can find it approximately by 
analyzing the Fourier transform of a typical local magne- 
tization profile of the rare region. Consider a rare region 
of linear size L (located at the origin) . Following Millis et 
al. 24 , the order parameter is approximately constant for 
r < L, while for large r > L, it decays as e~ r '^/r, where 
£ is the bulk correlation length. Taking the Fourier trans- 
form we find that F e (q) depends on e via the combination 
|Q — q| 3 log(e _1 ) only, where Q is the ordering wave vec- 
tor. Correspondingly, from Eq.(8), we find that the rare 
region contribution to the zero-temperature susceptibility 
in the quantum Griffiths phase can be expressed as 

Im X (q, u) cx M^sgnH X[(q - Q) 3 log^ 1 )] , (9) 

where X is a scaling function. The precise form of the 
logarithmic correction is difficult to find and beyond the 
scope of this paper. For random quantum Ising models, 
the susceptibility has the same structure as Eq.(9) [TT] . 
It is clear that the scaling function X will give only loga- 
rithmic corrections to the temperature dependence of the 
resistivity p in our further calculations. 

To minimize the resistivity functional (4), we need to 
make an ansatz for the distribution Close to an antifer- 
romagnetic quantum phase transition, the magnetic scat- 
tering is highly anisotropic because %(q, u>) peaks around 
the ordering wave vector Q. However, since we are inter- 
ested in a strongly disordered system, the low-temperature 
resistivity will be dominated by the elastic impurity scat- 
tering which is isotropic and redistributes the electrons 
over the Fermi surface. Correspondingly, we can use the 
standard ansatz 



$k«n-k. 



(10) 



where n is a unit vector parallel to the electric filed. Note 
that any constant prefactor in "I^ is unimportant because 
it drops out in the resistivity functional (4) and in the 
corresponding thermal resistivity functional (13). Then, 



after applying standard techniques 22 the magnetic part 



of the resistivity given in Eq. (4) becomes 



ApcxT 



d 3 < 



(n-q) 



~dw^Im X (q )W ). (11) 



Inserting the susceptibility (9) yields the rare-region 
contribution to the resistivity in the antiferromagnetic 
quantum Griffiths phase as 



Ap cx T 



(12) 



fol- 



Thus, the temperature-dependence of the resistivity „., 
lows a non-universal power-law governed by the Griffiths 
exponent A. 

Other transport properties. In the same way, 
we study other transport properties such as the thermal 
resistivity, the thermopower, and the Peltier coefficient. 
The variational principle for the thermal resistivity has 



the form 22 



VF[$k] = min 



J/($ k -%) 2 rk'dkdk' 



L (/v k ( £ 



(13) 



where \i is the chemical potential of the s-electrons. As 
long as impurity scattering dominates, we can use the 
standard ansatz for the variational function, 



<&k cx (gjj - p)n- k. 



(14) 



Then, following the calculation for the thermal resistivity 
outlined in Ref. |22| we obtain 



AW 



oc 



X 



1 

oj 2 



d 3 c 

v9 



(15) 



6k F 2 , 



2 

tt q 
3kp 2 



■T 2 



where kp is Fermi momentum of the s-electrons R In- 
serting the susceptibility (9) into (15), the temperature 
dependence of the thermal resistivity due to the spin- 
fluctuations in the Griffiths phase from the above equation 
is given by 



AW cx T 



A-l 



(16) 



The existence of an electric field E in a metal subject 
to a thermal gradient VT is called Secbeck effect and is 
characterized by the thermopower S which is defined via 
E = S VT. To calculate the thermopower, we analyze the 
Boltzmann equation (2) in the presence of both E and VT 
using the trial function 



cx ijin • k + ^(^k _ A*) n ■ k . 



(17) 



where 771 and 772 are variational parameters. Elastic im- 
purity scattering leads to the usual linear temperature de- 
pendence Simp T while the contribution due to the mag- 
netic scattering by the rare regions in the Griffiths phase 
reads 



AS" cx T 



A+l 



(18) 



2 Here, we have averaged over all directions of the vector n; this 
is sufficient to get the temperature dependence. 
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Fig. 1: (Color online) Schematic temperature-control parame- 
ter phase diagram of an itinerant antiferromagnet close to the 
quantum critical point. Our results apply in the Griffiths phase 
at low temperatures. 



Another transport coefficient called the Peltier coeffi- 
cient II characterizes the flow of a thermal current in a 
metal in the absence of a thermal gradient. It is related 
to the thermopower by II = S T. Correspondingly, the 
rare-region contribution to the Peltier coefficient has the 
form 



Alloc T x+2 . 



Discussion and conclusions. In summary, we 
have investigated the transport properties in the quan- 
tum Griffiths phase close to an antiferromagnetic quan- 
tum phase transition in a metallic system (see Fig. I). 
The rare-region contributions to electrical resistivity, ther- 
mal resistivity, thermopower, and the Peltier coefficient 
are characterized by non-universal power-laws in T which 
are controlled by the Griffiths exponent A. 

Our results have been obtained using the semi-classical 
Boltzmann equation approach. This approach is valid in 
Griffiths phase in which the system consists of a few lo- 
cally ordered rare regions in a non-magnetic bulk where 
the quasiparticles are well-defined. Sufficiently close to 
the actual quantum critical point (which is of infinite- 
randomness type) the quasiparticlc description may break 
down, invalidating our results. A detailed analysis of this 
question hinges on the fate of the fermionic degrees of free- 
dom at the infinite-randomness quantum critical point. 
This difficult problem remains a task for the future. 

We have used the standard isotropic ansatz (10,14) for 
the deviation of the electron distribution from equilibrium. 
This is justified as long as the rare-region part Ap(T) of 
the resistivity is small compared to the impurity part po- 
When A/9 becomes larger, the anisotropy of the scatter- 
ing needs to be taken into account. This can be done by 
adapting the methods of Rosch 25 to the situation at 
hand. 

We emphasize that our results have been derived for 
antiferromagnetic quantum Griffiths phases and may not 



be valid for ferromagnetic systems. The problem is that 
a complete theory of the ferromagnetic quantum Griffiths 
phase in a metal does not exist. In particular, the dynam- 
ical susceptibility still is not known. Correspondingly, the 
transport properties in ferromagnetic quantum Griffiths 
phases remain an open problem. 

Non-universal power-laws in a variety of observables in- 
cluding transport properties can also arise from a different 
physical mechanism far away from the magnetic quantum 
phase transition. In Kondo-disordered systems, the ex- 
istence of a wide distribution of local single-ion Kondo 
temperatures is assumed, this leads to the power-law sin- 



gularities 26 27 . This model was used to explain experi- 



mental results in some heavy fermion compounds such as 
UCu 4 Pd and UCu 3 . 5 Pdi. 5 



28 29 



Let us now turn to experiment. Unfortunately and 
somewhat ironically, all clear-cut experimental observa- 
tions of quantum Griffiths phases are in itinerant fer- 
romagnets rather than in antiferromagnets. However, 
quantum Griffiths effects have been discussed in the con- 
text of the antiferromagnetic quantum phase transition in 
heavy fermion systems (30 31 . One of the most striking 



predictions following from our theory is that the expo- 
nent characterizing the electrical resistivity should be less 
than one sufficiently close to the quantum phase transi- 
tion. There are several antiferromagnetic systems such 



(19) as CeC01.2Cuo.8Ge2 and Ce(Ru . 6 Rho. 4)2812 30 31 that 



show unusual power-law behaviour of the electrical resis- 
tivity with an exponent less than unity. The first system's 
resistivity increases with decreasing temperature. This 
is incompatible with our prediction and described by the 
Kondo model. The resistivity of the second compound 
decreases with decreasing temperature in agreement with 
our prediction. However, it is not clear whether this be- 
haviour is indeed caused by the quantum Griffiths phase. 
To establish this, one should measure various thermody- 
namics quantities as well as the transport properties and 
relate their low-temperature behavior. 
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